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Abstract 

The multirank separable kernels of the neutron-proton interaction for uncou- 
pled S and P partial waves (with the total angular momentum J =0,1) are 
proposed. Two different methods of a relativistic generalization of initially 
nonrelativistic form factors parametrizing the kernel are considered. Using 
the constructed kernels the experimental data for phase shifts in the elastic 
neutron-proton scattering for the laboratory energy up to 3 GeV and low- 
energy parameters are described. The comparison of our results with other 
model calculations are presented. 
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1. Introduction 

The problem of an adequate description of nuclear interactions arose 
many years ago. An important role in this task is played by the construc- 
tion of the nucleon-nucleon interaction. The simplest way to investigate such 
interaction is to describe properties of the elastic neutron-proton (np) scat- 
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tering and their bound state - deuteron. The latter can be considered through 
some reactions, such as the photodisintegration, the electrodisintegration etc. 

There are a lot of works devoted to the description of the deuteron. The 
first models were based on the nonrelativistic Shrodinger equation (see, for 
example, [H, 0]). In numerous papers mesonic exchange currents, relativistic 
corrections were investigated and then were added to nonrelativistic solu- 
tions These approaches could describe properties of the deuteron 
such as binding energy, magnetic and quadruple momenta, electromagnetic 
form factors, tensor polarization and so on. However, with increasing of the 
precision of experimental data and obtaining new data at larger energies it 
became evident that it was necessary to take into account relativistic effects 
more carefully. And the consistency of the consideration of the deuteron 
breakup reactions demands also the final state interaction (FSI) between the 
outgoing nucleons to be taken into account. 

One of the most consistent approaches is based on the solution of the 
Bethe-Salpeter (BS) equation [17|. In this case, we have to deal with a 
nontrivial integral equation. In addition to introduce the FSI of the final 
np pair the BS equation for the continuous state should be solved. There is 
no method to get its exact solution. So various approximations were worked 
out. At present, the most known and developed approaches are the so-called 



quasipotential 18]- 29 



They consist in the simplification of the equations under consideration 
by some assumption. In most cases it means the deliverance from the rel- 
ative energy, according to some ph ysic al reasons. For example, one of the 
particles is supposed on- mass-shell 27L or the time coordinates of particles 



are made equal [26j], etc. There is another approach called the light front 
(LF) dynamics which was successfully developed and applied to explain the 
electromagnetic properties of the deuteron. In this approach, the state vector 
describing the system is expanded by Fock components defined on a hyper- 
sphere in the four-dimensional space-time. The LF dynamics approach is 
intuitively appealing, since it is formally close to the nonrelativistic descrip- 
tion in terms of the Hamiltonian and state vectors can be directly interpreted 
as wave functions (see, for example, review 16j). The equivalence between 
LF dynamics and BS approaches was a subject of discussions presented in 



16l . |30j and references therein. 

An alternative approach based on the exact solution of the BS equation 
is to use the separable ansatz for the interaction kernel in the BS equation 
"3o| . In this case we can transform an integral equation to a system of linear 
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equations. Parameters of the kernel are fitted by the description of phase 
shifts for respective partial states and low-energy parameters. 

First separable parametrizations were worked out within nonrelativistic 
models. The form factors in the interaction kernel had no poles on a real axis 



in the relative energy complex plane |3ll . |32| . However, after the construction 



of a relativistic generalization such poles appeared [30( • In some cases they do 
not prevent to perform the calculations. As an example, such parametriza- 
tions are successfully used in the consideration of the deuteron photodisinte- 



gration 33], elastic electron-deuteron scattering, deuteron electrodisintegra- 
tion on the threshold, and the deep inelastic scattering 3(| . However, at high 
energies, one would have to deal with several thresholds corresponding to the 
production of one, two and more mesons of different types. This is clearly not 
feasible. A more practical approach is to employ phenomenological covariant 
separable kernel, which do not exhibit the meson-production thresholds, and 
can even be constructed in a singularity-free fashion, with the form factors 
chosen in the present paper and our Wick-rotation prescription. Thus, an 
accurate description of on-shell nucleon-nucleon data is possible, up to quite 
high energies. One then hopes that the used separable interactions also have 
a reasonable off-shell behaviour, so that realistic applications to other reac- 
tions can be done. The parametrization like that was proposed in [341] . In 



our previous works |35|, [36|, we constructed the one-rank interaction kernel 
of the same type. Due to the simplicity it works till the laboratory energy 
^Lab ~1 GeV except the simplest partial wave 1 P 1 + where it gives satisfactory 
results for all available experimental data. In the present work, we develop 
this approach increasing the rank of the kernels and trying to describe the 
data for the phase shifts of uncoupled partial states in the energy range up 
to 3 GeV taken from the SAID program (http://gwdac.phys.gwu.edu). In 
future we plan to use these kernels in relativistic calculations of the deuteron 
electro disintegration far from the threshold. 

It should be emphasized that in this paper the new kernel is fitted to de- 
scribe the neutron-proton elastic scattering data only. The description of the 
proton-proton scattering is a separate problem which requires specific meth- 
ods of calculation. Our aim was to construct the separable kernels suitable for 
consideration of the scattered 1 P\, 3 Po, 3 Pi, 1 So states of the neutron-proton 
system. 

The paper is organized as follows. In Section^ the general Bethe-Salpeter 
formalism is considered. The used separable kernel is described in Section 
El Section H] is devoted to the methods of a relativistic generalization of 
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nonrelativistic Yamaguchi- and Tabakin-type form factors. In Section [5J 
the pole structure of the obtained relativistic expressions is analyzed. The 
parametrizations for definite partial channels are presented in Section [61 In 
Section 0, the scheme of performing numerical calculations is offered. The 
obtained results and the comparison with other model calculations are dis- 
cussed in Section [SJ In conclusion, in Section EH the fields of application of 
the constructed kernels is briefly outlined. 

2. Bethe-Salpeter formalism 

Within the relativistic field theory, the elastic NN scattering can be de- 
scribed by the scattering T matrix which satisfies the inhomogeneous BS 
equation. In momentum space, the BS equation for the T matrix can be (in 
terms of the relative four-momenta p' and p and the total four-momentum 
P) represented as: 

T(p',p;P) = V(p',p;P) + ^ J d A kV(p',k-P)S 2 (k-P)T(k,p-P), (1) 

where V(p',p; P) is the interaction kernel and S^A;; P) is the free two-particle 
Green function 

S^\k-P) = (\P ■ 1 + k ■ 1 - m) {1) (\P ■ 1 - k ■ 1 - m) {2 \ 

7 are the Dirac gamma-matrices. The square of the total momenta s = 
ipx + P2) 2 and the relative momentum p = (pi — P2V2 [p' = (p[ — pQ/Z] are 
defined via the nucleon momenta pi, pi \p[, p' 2 } of initial [final] nucleons. 

To perform the partial- wave decomposition of the BS equation (JTJ), we 
introduce relativistic two- nucleon basis states \aM) = | vr, 2S+1 LjM), where 
S denotes the total spin, L is the orbital angular momentum, and J is the to- 
tal angular momentum with the projection M; relativistic quantum numbers 
p and 77 refer to the relative-energy and spatial parity with respect to the 
change of sign of the relative energy and spatial vector, respectively. Then 
the partial-wave decomposition of the T matrix in the center-of-mass frame 
(cm.) has the following form: 

T al 3,js(p',P;P(0)) 

= J2 (yaM(-p')U c ) a p ® (U c yl M (p)) 5j T ab (p' , |p'|;p , |p|; s), (2) 

JMab 
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where Uc = 27 2 7° is the charge conjugation matrix; the total momentum of 
the colliding nucleons in cm. is denoted by P( ). Greek letters (a, (3) and 
(7,5) in Eq.(T5]) refer to spinor indices and label particles in the initial and 
final states, respectively. It is convenient to represent the two-particle states 
in terms of matrices. To this end, the Dirac spinors of the second nucleon are 
transposed. At this stage T is 16 x 16 matrix in spinor space which, sand- 
wiched between Dirac spinors and traced, yields the corresponding transition 
matrix elements between SXJ-states. 

The spin-angular momentum functions y a M{p) are expressed in terms 
of the positive- and negative-energy Dirac spinors u^ 1 ^' 1 ^ 2 , the spherical 
harmonics Yhm L and Clebsch-Gordan coefficients Cj 1 mij 2 m2 : 

y,JM:LSp{v)Uc 

_ -L \^ (-tSpp r JM r Sm s 

^ / ^ 1 1 ^ T m r n ^ 1 1 

m L m s m 1 m2Pip2 

^(P)< (1) (P)< (2)T (-P). (3) 

The superscripts in Eq.flH]) refer to particles (1) and (2). To derive the matrix 
elements between a-states, the ortonormalization condition for the functions 
y a M{p') should be used: 

J dp P d(cos9 p ) Tr {^l M (p)^a'M'(p)} 

= Jdcp p d(cOs9 p )(yl M (p))p a (y a >M>{p))a0 = 8aa>5MM', (4) 

where partial states a and a' belong to the same partial channel. 

The partial-wave decomposition for the interaction kernel V of the BS 
equation (TjQ) can be written analogously to Eq.(T5]): 

= ^2(y a M(-p')U C )aP ® (Ucyt M (p))8y Vab{p' , \p'\]PQ, \p\] s). (5) 

abM 

Applying the condition (j3j), we can obtain a system of linear integral equa- 
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tions for the off-shell partial-wave amplitudes: 

Tabip'o, Ip'Upo, |p|;s) = Vah(p' , \p'\;p , |p|;s) (6) 



+00 

1 

471-3 

cd 



^ / <^o / k 2 <i|k| Kc(Po) Ip'I^o, |k|; s) 

„ ,1 J J 



x Scd(k , |k|; s) T dfe (fc , |k|;p , |p|; s), 

where the two-particle propagator S a b depends only on p-spin indices. 

We use the normalization condition for the T matrix in the on-mass-shell 
form for the singlet case: 

T a (s) =T a (0,p;0,p;s) = exp{zfr} sin 5 h (7) 

y'sy/s — Am z 

where p = |p| = — m 2 = \J mT^t, / 2 is the on-mass-shell momentum, 

m is a nucleon mass. In Eq.(j7]) I denotes s Lj states for simplicity. Low- 
energy parameters, the scattering length do, and the effective range r are 
derived from the expansion of the T- matrix into a series of p-terms, according 
to 



37 



p cot 6i(s) = -\ + ^p 2 + 0(p 3 ). (8) 

2 



To solve the equations for the T matrix and BS amplitude, we should use 
some assumption for the interaction kernel. 



3. A separable kernel 

We assume that the interaction kernel V conserves parity, total angular 
momentum J and its projection, and isotopic spin. Due to the tensor nuclear 
force, the orbital angular momentum L is not conserved. Moreover, the 
negative-energy two-nucleon states are switched off, which leads to the total 
spin S conservation. The partial-wave-decomposed BS equation is therefore 
reduced to the following form: 

Ti>i(p'o, Ip'I;Po, |p|; s) = V n (p' , \p'\;p , |p|; s) (9) 

+00 00 

y [ dk f k 2 d{k{ Vi'i'iPo, \P'\-M \k\;s) T^kp, |k| p , |p|; s) 
pi J M - + ie) 2 - k 2 



47T 3 
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where I = I' = J for spin-singlet and uncoupled spin-triplet states. 

Supposing the separable (rank N) ansatz for the kernel of the NN inter- 
action: 

N 

Vi'i(p' , \p'\;po, |p|;s) = x ij( s )9 [ i\p'o^ \p'\)gf(Po, IpI), (10) 

where the form factors gf represent the model functions, we can obtain the 
solution of equation fflOl) in a similar separable form for the T matrix: 

N 

Tin(p' , |p'|;pd, IpI; s ) = %( S )#F Vo> \v>'\)gf(Po, IpI), (h) 

where 

r,,( S ) = l/(A i r. 1 ( S ) + ^.( S )), (12) 

^' (s) = -^z2J dko J k rf|k| (TiT^^W^I' (13) 

Ajj(s) is a matrix of model parameters. 

The form factors used in the separable representation of the interaction 
kernel (fTUj) are obtanied by a relativistic generalization of the initially non- 
relativistic Yamaguchi-type functions depending on the three-dimensional 
squared momentum \p\. There are two methods to derive covariant relativis- 
tic generalizations of nonrelativistic form factors. They are considered in the 
next section. 

Calculating the T matrix we can connect the parameters of the internal 
kernel with observables. 



4. Methods of a covariant relativistic generalization 

In this section, two methods of a covariant relativistic generalization of 
the Yamaguchi- and Tabakin-type functions are presented. 

1. One of the common methods is to replace three-momentum squared by 
four-momentum squared: 

P 2 = ~pI + P 2 - (14) 

This formal procedure converts three-dimensional functions to covari- 
ant four- dimensional ones. 
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2. The other method is based on the introduction of the formal four- vector 
Q via the relative p and total P four-momenta of the two-body system 
by the following relation: 

Q = P-^P, (15) 

s 

with the total momentum squared s = P 2 . 

Note that in the two-particle center-of-mass system where P = (y/s, 0) the 
four- vector Q is defined by the components Q = (0, p) and, thus, 

p 2 = -Q 2 (16) 

can be formally converted to the Lorentz invariant. 

Let us consider the methods described above as applied to the nonrela- 
tivistic Yamaguchi-type function 

9 (|p|) = ^- (") 

In the first case, using the substitution ( Tl4l) we obtain the covariant func- 
tion in the form: 

9 P (P, p ) = 2 , „ 2 ^ — fl2 , ■ ■ (18) 

— p 2 + p 2 ~Pq + p 2 + p 2 + te 

In the second case, we use relation ffl6l) and obtain the function: 

^' P) = ^g^^p^- (19) 

The presented functions have rather different properties in the relative 
energy p complex plane in cm. The function g p has two poles on the real 
axis for p at ±a/p 2 + f3 2 =F «e while the function gq has no poles on it. 

In practical calculations of the reactions with the high momentum transfer 
the po integration can lead to singular expressions in functions g Pt Q on |p| 
or cos^p. This problem can be easily solved by calculating the |p| or cos# p 
principal value integral. However, another form of functions g Pt Q with odd 
powers in the denominator leads to nonintegrable singularities. Therefore, 
we introduce functions g Pt Q of type without poles on the real axis in the 
relative energy po complex plane. As an example of a function like that we 
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introduce the covariant form factors in the following form (see also section 3 
of H): 



a( v P) = Pc ~ p2 ±5> Pc-Pl + P 2 m 

9p^ r > ( p 2_ ^2)2 + ^4 (p2- p 2_ ^2)2 + «4' ^) 

and in the second case we use relation (fT6l) and obtain the function 

( p\~ Pc ~ Q 2 c - m - Pc + p2 foil 

^ J (Q2_ £2)2 + ^4 — * ( p2 + ^2)2 + ^4- ^ > 

We note that the function gq still has no poles on the po real axis while g p 
has poles at p : ±a/p 2 + (3 2 + ia 2 , ±a/p 2 + /3 2 — 2a 2 . 

The two methods of a covariant relativistic generalization described above 
can be investigated by solving the Bethe-Salpeter equation for specific partial 
states. 



5. Pole structure of the BS solution 

The solution of the Bethe-Salpeter equation with the separable kernel of 
interaction contains the function hij ( TTBl) . To simplify the investigation of 
the pole structure, let us consider for the one-rank kernel (i = j — 1) for 
single I state. Then the value to be calculated is h(s): 

(22) 

To obtain the function h(s), the two-dimensional integral on p and |p| 
should be calculated. To perform the integration over p , the Cauchy theorem 
is used. As it can be seen from Eq.()2"2"]). there are two types of singularities 
on the real axis in the po complex plane: one is poles of the function S(p; s): 

rf' 2) = ±v^/2 T E p ± ze, (23) 

and the other is poles of the function g(p). 
The function g p has four poles: 



Po 



(3,4) 



Po 



(5,6) 



ivV + Z^ + ia 2 , 

= ±y/p* + P 2 -ia 2 , (24) 
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and to perform the po integration, residues in three poles of Eqs. (l23[) and 
( I24p should be calculated. These calculations are performed analytically. 
This procedure is worthy of a special discussion. All poles and the contour 
of integration are pict ured in FigJTJ The idea how to choose the contour 



appeared owing to [38|, [39|. It consists in that the contour must envelope the 
poles from form factors which will be inside the standard contour after the 
a — > limit. "Standard" means the one used in the quantum field theory 
calculations with a propagator which has poles only on the real axis in the 
Po complex plane; one of them is rounded from below and the other, from 
above. So the path of integration is defined by an appropriate contour for 
the propagator. The calculation over the presented path leads to the pure 
real contribution from the form factor poles and, therefore, to the unitary T 
matrix. We also obtain a correct transition to ordinary form factors of type 
g ~ l/(pl — p 2 — I3 2 ) 2 in the a — ► limit. 

The function gq has no poles on the po real axis and, therefore, the only 
poles of Eq. (j251) should be taken into account. The result for h(s) can be 
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written as: 



2vr 2 J * ^^-2E p + ie 
This equation formally coincides with that could be obtained within the 



, , , ■ / 2 ji i 9q(®, IpI) 2 / oc; x 



Blankenbeckeler-Sugar-Logunov-Tavkhelidze (BSLT) approximation [181 . Il9 
which consists in replacing the Green function in Eq. (1221) by the expression 

Sbslt(p; s) = -2m{^s~ - 2E P + ze)" 1 ^). (26) 

Although the solutions of the equation with functions $q and within the 
BSLT approximation coincide in cm., the difference becomes evident when 
the reaction with the two-particle system is considered. In that case, the 
arguments of the function gQ are calculated with the help of the Lorentz 
transformations in the system different from cm. The functions become 
similar to g p but with a more complicated dependence of the p argument. 
In our calculations we prefer to use the first method of relativization. The 
reason is in that reaction we are planning to consider in future (namely, 
electrodisintegration) these transformations could lead to the appearance of 
additional poles in the calculated expressions. However, in this work the 
comparison of results for phase shifts and low-energy parameters is presented 
for both cases. 

In following two sections some separable presentations of the interaction 
kernel for the partial waves with J = 0, 1 are considered. Form factors are 
constructed by the relativization procedure, according to the first method. 
The functions with Q can be obtained from them by the change p 2 — » Q 2 . 

6. Separable presentations of the kernel 

We consider the partial states using the separable kernels with modified 
Yamaguchi-type functions of N rank (MYiV, MYQiV). 

6.1. Two-rank kernel for P states: 3 Pq, l P\, 3 P\ 

For the description of the P partial waves the two-rank separable kernel 
of interaction with the modified Yamaguchi-type functions is used. Its form 
factors are written as: 

s'PiP) = V f° + J* * • < 27 > 

(Po - P - Pi) 2 + ot\ 



9 [ P(P) 



V(-P 2 + P 2 ) 3 (Pc2 ~ Pi + P 2 

{{p 2 - p 2 - (3 2 ) 2 + aj) 2 
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6.2. Three-rank kernel for 1 Sq state 

The investigation we performed demonstrates a bad description of phase 
shifts for 1 Sq partial state by the two-rank kernel. It is not much better than 
by the one-rank kernel |35] . So for this special case the three-rank interaction 
kernel was elaborated. Its form factors have the following form: 



[SI, ^ _ (gel ~ V\ + P 2 
[S],s _ (P 2 Q-P 2 )(Pc2-P 2 q + P 2 ) 2 



((p 2 , - p 2 - & 2 ) 2 + «t) 2 



[S], x (go - P 2 ) 

#3 (P) " 



4 ' 



(p 2 , - p 2 - & 2 ) 2 + «: 

The functions are numerated by angular momenta L = 0([S]), 1([P]). The 

[PI * c 

numerator with p c in ^ an d g[ 2 is introduced to compensate an additional 



dimension in the denominator to provide the total dimension as GeV 2 35 



7. Calculations and results 

Using the np scattering data we analyze the parameters of the separable 
kernels distinguishing three different cases: 

1. There are no sign change in phase shifts or bound state ( 1 -P 1 f , 3 P^~ 
partial states). In this case 

Ajj(s) = Xij = const. (29) 

This is sufficient for most of the higher partial waves. 

2. One sign change and no bound state ( 1 S' f , 3 P t ~ partial states). In this 
case the energy- dependent expression for Xi(s) is used (see jiij and 
references therein): 

Kj(s) = (s Q - s)Xij, (30) 

Here the parameter so is introduced to reproduce the sign change in the 
phase shifts at the position of the experimental value for the kinetic 
energy TLab where they are equal to zero. It is added to the other 
parameters of the kernel. 
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Figure 2: Contour for integration over po after the Wick rotation. 



The calculation of the parameters is performed by using Eqs.fllJ), (jHJ) and 
expressions given in two previous sections to reproduce experimental values 
for all available data from the SAID program ( |http: //gwdac.phys.gwu.edu[ ) 
for the phase shifts. The the low-energy scattering parameters are taken 
from 40| . 

The calculations are performed in two independent ways. The first one 
is based on using the Cauchy theorem for the integration over p compo- 
nent of the nucleon four-momentum. The integration over |p| is performed 
numerically. The second one is based on using the Wick rotation 39]. In 
this case the contour of integration Fig{T] is deformed as depicted on Figj2j 
All integrals are calculated numerically with the technique elaborated in the 
paper 
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Now we find the introduced parameters of the kernel: 

1. For P waves the minimization procedure for the function 



X 2 = 5> eXP ( s *) - ^))7(A5 CXP (^)) 2 (31) 
i=i 
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is used. Here n is a number of available experimental points. 
2. For the l S$ wave the values of the scattering length a is also included 
into the minimization procedure 

n 

X 2 = $> GXP ^) - ^))7(A^ CXP (^)) 2 + K xp - a) 2 /(Aa cx P) 2 . (32) 

i=l 

The effective range r$ is calculated via the obtained parameters and compared 
with the experimental value rg Xp . 

The description of P waves by the two-rank kernel is denoted by MY2 for 
the first case of a relativization procedure; MYQ2, for the second one. For 
the 1 Sq partial state the notation MY3 and MYQ3, respectively, are used. 

The calculated parameters of the considered kernels are listed in Tables [1] 
and [2] (here the values of s are presented, too). In Table [HJ the calculated 
low-energy scattering parameters for the 1 Sq wave are compared with their 
experimental values. 

In FigsOHHl the results of the phase shift calculations are compared with 
experimental data (the used notation is described in the following Section 



IE]) and two alternative descriptions by CD-Bonn [42| and SP07 43j. In the 
discussion of the 1 Sq channel the nonrelativistic Graz II 31], as an alternative 
model with a separable kernel, is also presented. 



8. Discussion 

In this section, the review of the results of our calculations with two 
methods of relativization is performed. 

In. FigE] we can see that all of the calculations including nonrelativistic 
CD-Bonn give a quite good description of the 1 P 1 + state. In our previous 



works [35|, [36[, it was shown that the one-rank kernel is already sufficient 
to reproduce the phase shifts in this case. The reason is the simplicity of 
their behavior and the narrowness of the energy interval where they are 
known. So if the simplicity of the model to perform calculations in the 1 P 1 + 
state is preferable, it is sufficient to use the one-rank kernel. For the energies 
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Table 1: Parameters of the two-rank kernel with modified Yamaguchi functions for P 
waves. 





MY2 






3p+ 
^0 




An (GeV 4 ) 
A12 (GeV 4 ) 
A 22 (GeV 4 ) 


0.05412952 


-923.8881 


0.06125619 


1.925 


-102.0961 


2.068215 


7.975 


4.346553 


24.48148 


Pi (GeV) 


0.1244769 


0.958602 


0.1224502 


fa (GeV) 


0.6228701 


1.897255 


0.5822389 


ai (GeV) 


0.2 


0.759970 


0.2107709 


a 2 (GeV) 
Pc2 (GeV 2 ) 
s (GeV 2 ) 


0.5984991 


0.687087 


0.5927882 


1.035 


-133.2385 


0.9476951 




3.8682 






MYQ2 






3 p+ 
M) 


3p+ 


An (GeV 4 ) 
A12 (GeV 4 ) 
A 22 (GeV 4 ) 


0.3086486 


-55.88270 


0.07648584 


1.606382 


-763.1649 


1.567463 


-5.797411 


-5325.327 


21.25497 


Pi (GeV) 


0.2150637 


0.5008744 


0.2109235 


Pi (GeV) 


0.9582849 


2.4269161 


0.4260685 


ai (GeV) 


0.2 


0.6863803 


0.2 


a 2 (GeV) 
Pc2 (GeV 2 ) 


0.2 


0.2 


0.2 


12.47736 


5.866078 


0.007380749 


so (GeV 2 ) 




3.8682 
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Table 2: Parameters of the three-rank kernel with modified Yamaguchi functions for 1 

state. 

MY3 



An (GeV 2 ) 
Ai 2 (GeV 2 ) 
A13 (GeV 2 ) 
A 22 (GeV 2 ) 
A 23 (GeV 2 ) 
A33 (GeV 2 ) 
so (GeV 2 ) 


-0.2922173 


ft (GeV) 


0.7018063 


-3.953624 


ft (GeV) 


4.381178 


1.035416 


ft (GeV) 


1.137604 


-11268.52 


ai (GeV) 


1.297282 


331.1130 


a 2 (GeV) 


4.612956 


70.37369 


a 3 (GeV) 
p c i (GeV 2 ) 
p C 2 (GeV 2 ) 


0.6752485 


4.0279 


38.20462 






34.53211 



MYQ3 



An (GeV 2 ) 
A12 (GeV 2 ) 
A13 (GeV 2 ) 
A 22 (GeV 2 ) 
A 23 (GeV 2 ) 
A33 (GeV 2 ) 
so (GeV 2 ) 



723.7399 
550.4827 
-1583.031 
132.7836 
-14.26609 
-26005.23 
4.0279 



ft (GeV) 
ft (GeV) 
ft (GeV) 
en (GeV) 
a 2 (GeV) 
a 3 (GeV) 
Pa (GeV 2 ) 
Pc2 (GeV 2 ) 



2.463736 
1.266692 
7.714815 
3.752405 
2.117966 
1.338336 
-44.84322 
247.2476 



Table 3: The low-energy scattering parameters and for the singlet (s) x Sq wave. 





a s (fm) 


r 0s (fm) 


MY3 


-23.750 


2.70 


MYQ3 


-23.754 


2.78 


Experiment 


-23.748(10) 


2.75(5) 
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^Lab > 11 GeV the resulting functions become very different. To make choice 
in favour of one of the parametrizations, it is necessary to have experimental 
data for larger energies. 

In FigJH the results of the calculations for the 3 P + partial state are de- 
picted. The comparison of our and other model calculations demonstrates a 
reasonable agreement with the experimental data in the whole range of en- 
ergies except the nonrelativistic CD-Bonn which works till T^b ~ 0.75 GeV. 
The good description of the phase shifts in 3 P + can also be archived by using 
the one-rank interaction kernel with extended Yamaguchi-type form factors 



sec 



35l . |36j). As in the previous case, the choice of the kernel is defined, by 



reasons of convenience, as applied to a specific problem. 

The phase shift calculations for the 3 P 1 + partial state are presented in 
Figj5j All models except CD-Bonn give various but acceptable within the 
limits of error results. The CD-Bonn works for TL a b ^ 0.66 GeV. As for our 
previous one-rank kernel, the increase of the rank allows us to improve the 
description and embrace all experimental data. 

Our general conclusion about the description of P states is that using the 
two-rank kernel allows to reproduce the phase shifts in good agreement with 
experimental data. If the simplicity is the main requirement of performing 
calculations, the one-rank kernel for 1 P 1 + and 3 P + states and two-rank one 
for 3 P X + can be used. Here we talk about phases in the whole energy range 
for TLab- However, if the use of a common interaction kernel is preferable, 
then calculations with two-rank interaction kernel should be done. 

From FigJHl where the results of calculations for 1 Sq are presented, it 
can be seen that CD-Bonn is quite good till T^b ~ 0.66 GeV. All the other 
results are in agreement with measured phase shifts except Graz II which 
works for T Lab ^ 2 GeV. We succeed in an acceptable description at the cost 
of increasing the kernel's rank. Kernels of lower ranks are proper only for 
energies TL a b ^1 GeV. 

From the presented figures the similarity of the calculations with the MY 
and MYQ form factors can be noted. Thus, as for the description of phase 
shifts and low-energy parameters there is no difference which variant of a 
relativistic generalization to choose. The choice should be dictated by the 
convenience of performing calculations within some concrete problem. 
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Figure 6: Phase shifts for the 1 Sq wave. 
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9. Conclusion 

Using the multirank kernels (two-rank for P waves, three-rank for the 
1 Sq partial state) we have constructed an adequate description of all ex- 
istent experimental data for phase shifts taken from SAID and low-energy 
parameters with capable accuracy. 

The results for two different methods of a relativistic generalization of 
initially nonrelativistic Yamaguchi-type form factors were considered. As 
it was shown they lead to slightly different descriptions of phases and low- 
energy parameters. Hence, the choice of the concrete form of functions for 
performing calculations of any process is governed only by convenience. 

In spite of the fact that the model functions have a simple form there are 
quite a few parameters in the description of the data. This is necessitated 
by introduction of an additional parameter a so that integrands containing 
form factors of the separable kernel could not have poles. In particular, using 
this type of kernels will make the numerical calculations of the electrodisin- 
tegration far from the threshold possible without resorting quasipotential or 
nonrelativistic models. 
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